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Abstract 

In the presence of a chemical potential, the physics of level crossings leads to singularities at zero tem- 
perature, even when the spatial volume is finite. These singularities are smoothed out at a finite temperature 
but leave behind non-trivial finite size effects which must be understood in order to extract thermodynamic 
quantities using Monte Carlo methods, particularly close to critical points. We illustrate some of these is- 
sues using the classical non-linear 0(2) sigma model with a coupling (3 and chemical potential /i on a 2 + 1 
dimensional Euclidean lattice. In the conventional formulation this model suffers from a sign problem at 
non-zero chemical potential and hence cannot be studied with the Wolff cluster algorithm. However, when 
formulated in terms of world-line of particles, the sign problem is absent and the model can be studied 
efficiently with the "worm algorithm". Using this method we study the finite size effects that arise due to 
the chemical potential and develop an effective quantum mechanical approach to capture the effects. As a 
side result we obtain energy levels of up to four particles as a function of the box size and uncover a part of 
the phase diagram in the fi) plane. 

PACS numbers: ll.lO.-z, 02.70.Ss, 05.30.-d, 75.10.Hk 



*Electronic address: jdebasish @ theory, tifr.re s .in, sch @ phy. duke .edu 



1 



I. INTRODUCTION 



Understanding the phase diagram of quantum chromodynamics (QCD) as a function of tem- 
perature T and baryon chemical potential /i is an active area of research. Although much is known 
about the physics at = from lattice QCD calculations jl, 2], there is controversy of what 
might occur at non-zero jj, and small values of T ||3|-l5|]. Due to the sign problem, which arises in 
all current formulations of lattice QCD at non-zero /x, it is impossible to perform first principles 
calculations to settle the controversy today. Most of our knowledge of the (T, ji) phase diagram 
of QCD is based on models that are motivated from universality and solved using mean field the- 
ory. Can at least some of these models be studied from first principles? For example, recently a 
Landau-Ginzburg approach was used to uncover parts of the phase diagram of QCD where the low 
energy physics is described by bosonic excitations [|6Q. In these regions it should be possible to 
construct bosonic effective field theory models that share the same symmetries, low energy physics 
and possibly the phase transitions as QCD. It would be interesting to study these models from first 
principles. Unfortunately, sign problems also arise in bosonic field theories in the presence of a 
chemical potential when formulated in the conventional approach. For this reason not many first 
principles studies of field theories with a chemical potential exist. However, many of these sign 
problems are solvable today and thus allow us to explore the physics of a chemical potential from 
first principles. It may be useful to study these simpler field theories before attempting to study 
QCD. 

One of the simplest examples of a relativistic bosonic field theory is the classical non-linear 
0(2) sigma model on a cubic lattice which has been studied extensively in the context of superfluid 
transitions using the efficient Wolff cluster algorithm [|7D. The phase transition is between two 
phases: an 0(2) symmetric phase and a phase where the symmetry is spontaneously broken. 
Close to the phase transition the low energy physics is described by an interacting quantum field 
theory of massive charged bosons in the symmetric phase and of massless Goldstone bosons in 
the broken phase. At the critical point the low energy physics is scale invariant and the critical 
behavior belongs to the three dimensional XY universality class. 

Since the model contains an exact 0(2) global symmetry, one can also introduce a chemical 
potential /i that couples to the corresponding conserved charge. This chemical potential helps one 
study the "condensed matter" composed of the fundamental boson present in the theory. When fi ^ 
0, the action in the conventional formulation becomes complex and Monte Carlo algorithms suffer 
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from a sign problem exactly like in QCD. Not surprisingly, the phase diagram of the condensed 
matter arising in the 0(2) non-linear sigma model has not been studied from first principles. On the 
other hand non-relativistic bosonic lattice models, especially in the Hamiltonian formulation have 
been studied for many years by the condensed matter community. Here one naturally constructs 
the field theory with bosonic world lines and there is no sign problem when one introduces a 
chemical potential. Thus, it is natural that a world-line approach could solve the corresponding 
sign problem even for a relativistic field theory. This was shown explicitly for both the linear 
sigma model [8] and the 0(2) non-linear sigma model [Q]. 

While the world-line representation for bosonic lattice field theories was well known for many 
years, the main advance in the field that improved our ability to perform a first principles calcula- 
tion in the presence of a chemical potential, was the discovery of an efficient Monte Carlo algo- 
rithm called the "worm algorithm" lioll . Variants of this algorithm in the name of "directed loop 
algorithm" [llll . 11211 have been used to solve a variety of models that arise in the strong coupling 
limit of lattice gauge theories [13-1^. The worm algorithm has also been found to be an efficient 
approach to study a wider class of fermionic field theories in the loop representation in two dimen- 
sions where fermion sign problems are absent [|l7l- l20|l and weak coupling Abelian lattice gauge 
theory ll^ . A combination of the worm algorithm and the determinantal algorithm was recently 
developed to solve the lattice Thirring model in the fermion bag formulation in higher dimensions 
II22I1 . All these developments should allow us to explore the physics of a chemical potential using 
first principles in a variety of lattice models with interesting symmetries. 

In this work, we explore the 0(2) non-linear sigma model in the presence of a chemical poten- 
tial and show that interesting finite size effects naturally arise due to the level crossing phenomena. 
Understanding these effects is important to extract the thermodynamic limit and thus uncover the 
(/3, yu) phase diagram, where (3 is the coupling and yU is the chemical potential. Our studies should 
be useful for future studies since the finite size effects we uncover is a universal feature. Our 
work also provides accurate results that can be used to compare with results from other methods, 
like the complex Langevin method, which are being explored as a solution to the sign problem 
in general [12314250. Our work is organized as follows: In section HI] we discuss our model and 
observables in order to set up the notation. In section |lll] we discuss the finite size effects that arise 
in the presence of a chemical potential and develop an effective quantum mechanical description 
that captures these effects in section |IVl Sections |V] and |VI] contain our results obtained using the 
worm algorithm. In particular we show that the observed finite size effects are described well by 
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the effective quantum mechanical description of section |IVl In section IVIIl we discuss the phase 
diagram of the 0(2) model which emerges from our work. 



II. MODEL AND OBSERVABLES 

The action of the 0(2) non-linear sigma model on a lattice with a finite chemical potential that 
we study here is given by 

S" = -/3 ^ I COs{6a: - " i^^5a,t) | , (1) 

where x is the lattice site on a three dimensional cubic lattice, a = 1,2 represent the spatial 
directions and a = t represents the temporal direction. We will use L to represent the spatial size 
and Lt the temporal size and assume periodic boundary conditions. The constant /3 plays the role 
of the coupling. The chemical potential /i is introduced in the standard way and couples to the 
conserved charge of the global 0(2) symmetry [I26]. When /i 7^ the action becomes complex 
and Monte Carlo algorithms to generate configurations [9] that contribute to the partition function 

[de,]e-^, (2) 



suffer from a sign problem. In particular the Wolff cluster algorithm [|71] is no longer useful at 
non-zero chemical potential. Hence the phase diagram of the model in the (/3, n) plane remains 
unexplored. 

It is possible to avoid the sign problem if we rewrite the partition function in the world-line 
representation [i9J . Using the identity 

00 

exp{cos^}= 4(/3)e*'=^ (3) 
fc=— 00 

where Ik is the modified Bessel function of the first kind, on each bond (x, a), and performing the 
angular integration over Or^ the partition function can be rewritten as 

^ = E n {KMe'-'-'''-'-] S(^Y.(^x,a-K-a,a)), (4) 

[k] X a 

where the bond variables k^,^a describe "world-lines" or "current" of particles moving from lattice 
site X to the site x + a and take integer values. A configuration of these bond variables, denoted 
by [k], is thus a world-line configuration. The global f/(l) symmetry of the model is manifest 



in the local current conservation relation represented by the delta function. In other words any 
particle that comes into the site must leave the site due to current conservation. In this world-line 
formulation the partition function is a sum over explicitly positive terms even in the presence of 
/i. Details of the "worm algorithm" that we have developed to update the world-line configuration 
[k] is described in appendix lAl 

We focus on four observables in this work: 

1 . The average particle density p: 



xgtimeslice 

The average particle number is then given by (A^) = pL^. 
2. The particle number susceptibility k: 



Note thsitK = Lt/L^{N^). 
3. The superfluid density (or particle current susceptibility) ps'. 




^ a=l,2 \ X 

The superfluid density is known to be ps = l/Lt{W'^) where W is the spatial winding 



number of particles 12711 . We define (Ng) = LF'Ps as the number of particles that are in the 



superfluid phase in a finite system. 
4. The condensate susceptibility X- 

X = J](e^^^e-^^^) (8) 

y 

The first three observables are "diagonal" observables and can be measured on each world-line 
configuration and then averaged. The condensate susceptibility x o^i the other hand is a "non- 
diagonal" observable, but it can be related to the size of each worm update as discussed in ap- 
pendixlAl We discuss some tests of the algorithm in appendix|Bl In particular we have been able 
reproduce earlier results of the 0(2) non-linear si gma model at p = 0. One of these results is 



the estimate of the critical coupling /3c = 0.45421 [|28|1 . For jS > l3c the 0(2) symmetry is spon- 
taneously broken, while for (3 < (3^ the model is in the symmetric phase. Our tests also show 
agreement with exact calculations on a 2 x 2 lattice at non-zero p. 



III. FINITE SIZE EFFECTS 



A good understanding of finite size effects is important for extracting thermodynamic results 
from numerical computations. This is particularly true close to a second order critical point where 
correlation lengths diverge. While developing a theory of finite size effects, one usually assumes 
Lt = where z is called the dynamical critical exponent of the problem. Such a choice makes 
calculations natural. In a relativistic theory since 2; = 1, it is natural to choose Lt = L. Most 
studies of the critical behavior at /3 = /3c and /i = make use of this choice. On the other hand, in 
the presence of a chemical potential, since the low effective energy theory is non-relativistic, one 
expects z = 2 and Lt = is a more natural choice. However, in our work we have found that 
even with this choice the finite size effects close to the critical point are non-trivial in the presence 
of a chemical potential. In fact observables always show clear "wiggles" and cannot be fit to a 
simple power law that one expects near the critical point. In order to demonstrate this feature, in 
Fig. [T] we plot the behavior of the average particle density p as a function of fi for L = 12 and 
Lt = 144 at (3 = 0.43. From mean field theory we expect 



close to /i = in the thermodynamic limit. In Fig. \T\ we observe that p is indeed zero for 
/i < 0.27 and begins to increase for /i > 0.27. But the increase, although roughly linear close 
to /ic as expected, shows clear "wiggles" when 0.27 < /i < 0.38 and only for ji > 0.38 the 
"wiggles" disappear. The region between 0.27 < ji < 0.38 has been enlarged in the left inset 
in order to enhance the observed "wiggles". In the right inset we fix /i = 0.32 and plot p as a 
function of L assuming Lt = L'^. Again the data shows clear oscillations whose origin may seem 
a bit mysterious. In order to avoid these oscillations one has to go to much larger L at a fixed 
value of fi. However, since Lt scales like going to larger lattice sizes is more difficult than in a 
relativistic theory. For this reason, we believe it may be useful to develop a different type of finite 
size analysis. 

As we will argue below, the strange finite size behavior is the result of energy levels cross- 
ing each other due to the chemical potential. Since the particle number is a conserved quantum 
number, energy levels with different particle numbers can cross each other at critical values of the 
chemical potential. Similarly at a fixed chemical potential, the changes in spatial size can also 
cause these energy levels to cross each other. These level crossings lead to singularities at low 




(9) 
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FIG. 1: Plot of p as a function of /i for /3 = 0.43, L = 12 and Lt = 144. The data shows clear "wiggles" 
at small values of p which disappears for larger values. The left inset magnifies the region of the "wiggles". 
The right inset shows the plot of p as function of L at p = 0.32 assuming Lt = which also shows clear 
non-monotonic behavior. 

temperatures (large L^) in a finite spatial volume (fixed L). While these singularities are smoothed 
out at finite Lt, they leave behind non-trivial finite size effects observed above. While it may still 
be possible to develop practically useful finite size scaling relations using Lt = L^, we find it 
natural to consider a finite size scaling theory for quantities as a function of Lt and n for a fixed 
value of L close to the critical values of ji where energy levels cross each other. As we discuss 
below this leads to an effective quantum mechanics problem. The finite size effects studied here 
have been observed earlier in the context of quantum spin-systems in a magnetic field Jig], but 
they were not analyzed using the techniques we introduce below. 
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IV. EFFECTIVE QUANTUM MECHANICS 



At a fixed value of L for sufficiently large Lf, it must be possible to map the lattice field theory 
problem to an effective quantum mechanics problem where only a few low energy levels play an 
important role. Let us label these energy levels by |A^, k) and the energy eigenvalues by Ej^^\ 
where N = 0, 1, 2, ... represents the particle number sector of the energy level and k represents 
"other" quantum numbers. The levels and the energies depend on L and jS. The partition function 
of the problem may be written as 

k,N 

Using this effective quantum mechanical description we can in principle find the Lt and n depen- 
dence of various quantities. However, for the analysis to be practically useful we need to assume 
that only a few energy levels are important. If we assume that /i is close to a critical value Hc 
where level crossing phenomena occurs, then for large enough Lt one might expect the physics to 
be dominated by just two levels. In this approximation we will derive the Lt and dependence of 
all our observables. 

In a given particle number sector can assume Ef) < < ) < ... without loss of 

generality. However, in this work we will also assume that£;J°^ < E^^^ < E^^\.. which means that 
it always costs energy to add a particle into the system. While this is not necessary it is precisely 
the situation we encounter in this work and simplifies our analysis. With these assumptions it is 
easy to argue that close to the critical chemical potential where the particle number changes from 

to + 1 we can approximate the partition function to be 

Here we have assumed all higher energy states will be suppressed exponentially at large Lt. It 
is easy to verify that jii^^ = E^f^'^^^ — E^,^^ is the critical chemical potential where the average 
particle number changes from N to N + 1. Below we discuss the Lt and n dependence of each 
observable when /j, a; iJ,i^\ 
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A. Particle Number 



We first consider the average particle number (N). When A^^^ — 1^ — A*c^^ is small and Lt is 
large we can write 

We will demonstrate later that our data fits very well to this simple one parameter fit and we are 
able to extract jii^^ very accurately for all L < 20 for a variety of values of ^. 



B. Number Susceptibility 

Next we discuss the number susceptibility k = Lt/ LF'{N'^). We now obtain 

1 + e^M Lt 

The value of ixi'^ is the same as obtained from the average particle number. So this observable 
has no new free parameters. 

C. Superfluid Density 

In the effective quantum mechanical description the superfluid density is given by 

e-(^*-*)^OV*^0 I (14) 



= - / Tr ( f 
Z Jo \ 



where 

= -^^Ji(xi,X2) (15) 

is an operator in the Hilbert space made up of the conserved current operator Ji{xi,X2) in the 
direction i at the site with coordinates (xi, X2). Note the sum is over the surface perpendicular to 
the direction of the current. Since it is a conserved current it does not matter which surface one 
chooses. Now if we introduce a complete set of energy eigenstates we get 

^ fl_e-(4"''-4"^-("'-")M)i*'i 
Ps = ^ y e-(<'-"'')^* y \{n,k\0\n'k')\^^^, L (16) 

First we note that (n, k\0\n'k') oc (5„„/, since the current operator commutes with the particle 
number operator and hence does not change the particle number. Further, as before we assume 



only two low lying energy levels are important in the partition function when /i ^ /li^^ . Then the 
[A;, n] sum is replaced by A; = and n = A^, + 1. Hence we obtain 

where 

f 1 - e-^<'-0 

Po = EK^-Q|Q|^^^) r . ^ withA^f^^^^f^^-Ef). (18) 

and pi is obtained by replacing N with + 1 in the above expression. Note that the sum over 
k' does not contain the k' = sum because {N, 0|O|A^0) = since O is a current operator and 
the ground state is rotationally invariant. Thus, the expression for ps contains two new parameters 
since pc^^ has already been encountered before. 



D. Condensate Susceptibility 

The expression for the condensate susceptibility can also be obtained using eq.(fT6l) if the oper- 
ator O is replaced by 

= |5^cos(M- (19) 

The matrix element {n,k\0\n'k') is non-zero only when n' = n + 1 or n' = n — 1. When 
fi ^ p[^\ the [k, n] sum is again dominated by E^'^^ and E^f^^^K However, in the present case 
the p dependence also enters the k' sum. In the limit as A^!^^ — > and — )• oo the A;' = term is 
singular while the other terms are not. Separating the singular term from others we find 

X = 2\{N,0\O\N+l,0)\' ^ 



Ar(l + eAr^^.) 

(1 



+ .Vo (Ai^r^)-A(^)) 

^a^> S *"^WtW- '''' 

Since IAI/^-*! is assumed to be much smaller than all AE^f^ and AeI^'^^\ at large Lf the expo- 
nentials in the k' sum can be dropped. If the remaining terms are expanded in powers of A^^^ we 
find 

_ Xo(e<"^"- - 1)/Ar + (Xi + X2Ar + -) + (x; + X^.A^^) + ..■)e-^^"'^^ 
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We find that our data fits well to this expression truncated at the quadratic order in A/t ' , which 
means we have seven new parameters in our fit. However, most of these parameters are not de- 
termined reliably and contain large systematic errors. The only parameter that can be determined 
reliably is xo and this is what we quote as a result from our analysis. 



V. RESULTS 



We have performed extensive calculations at /3 = 0.43, 0.50 and 0.20. These values of (3 are 
chosen so that two of them are close to the critical coupling (3c = 0.45421 on either side and one 
is far from it in the massive (disordered) phase. In this section we present fits of our results to 
the effective quantum mechanics description discussed above. As mentioned earlier, the effective 
description becomes useful only in the limit of small temperatures where excitations to higher 
energy levels can be neglected. Since the spacing between energy levels decreases with increase 
in volume, our approach works best on small spatial volumes. However, thanks to the efficient 
worm algorithm, we have been able to extract parameters of the effective quantum mechanics up 
to L = 16 at /3 = 0.43 and /3 = 0.5. Although this lattice size is small compared to normal studies 
of bosonic lattice field theories, it still allows us to perform a useful study of the L dependence 
of the physics and draw quantitative conclusions about the thermodynamic limit. At (3 = 0.2 we 
observe that the energy levels are more densely packed and we are able to compute quantities only 
up to L = 8. 

We first consider L = 2 and vary Lt in the range 40 < < 200 which is easy due to the 
small lattice size. In Fig. [2] we plot the behavior of our four observables as a function of the 
chemical potential at Lt = 100 and (3 = 0.43. Note that the particle number increases in steps 
of one at critical values of fi. This means energy per particle of the ground state in every particle 
number sector increases with the number of particles. In other words the particles repel each other. 
Thus, small systems containing particles of the non-linear sigma model will show the phenomena 
similar to Coulomb Blockade observed in nanoscale systems [|30D. By fitting the data at L = 2 
and /3 = 0.43, 0.50 and 0.20 we have extracted the parameters yUc^^ po, pi and Xo for ^ = 0, 1, 2 
and 3. These are tabulated in Tab. HI In order to show the goodness of our fits, in Fig. [3] we plot 
the behavior of (N) and x for values of p close to the transition between the = and A^ = 1 
sector for different values of Lt at (3 = 0.43. The solid lines represent the fit functions using the 
parameter values given in table IH Note that all the computed values of (A^) shown in the left plot 
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FIG. 2: The four observables as a function of fi up to four particle excitation. The data shown is for L = 2, 
Lt = WO and (3 = 0.43. The solid Unes are fits to the effective quantum mechanical description. 

of Fig. [3]can be fit with just one parameter namely fxi^^ . 

We have repeated the above analysis at larger values of L. We find the physics remains qualita- 
tively similar to the L = 2 case. In particular the average particle number jumps by one at critical 
values of fi. In Fig. |4] we show the average particle number as a function of /i at different values 
of L at /3 = 0.43 and (3 = 0.50. In Figl5]we plot the average particle number at L = 8 for all 
three values of /3 at different values of Lt close to fii^^ for = 1,2, 3. The effective quantum 
mechanics description continues to fit all our data well as long as /i is close to the critical values 
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N 




Po 


Pi 


Xo 


xVdof 


/3 = 0.43 





0.82444(3) 


0.2114(2) 


0.2713(2) 


2.349(4) 


1.23 


1.43 


0.70 


0.58 


1 


1.22462(2) 


0.2715(2) 


0.2983(2) 


2.817(4) 


0.52 


0.50 


1.04 


0.81 


2 


1.51567(2) 


0.298(1) 


0.312(1) 


3.013(4) 


0.73 


0.70 


1.20 


1.28 


3 


1.74436(2) 


0.313(1) 


0.321(1) 


3.121(4) 


1.02 


1.02 


0.78 


1.14 


/3 = 0.50 





0.63275(3) 


0.288(1) 


0.336(1) 


2.63(1) 


0.58 


0.56 


1.33 


0.68 


1 


1.05865(2) 


0.339(1) 


0.362(1) 


2.957(9) 


0.72 


0.69 


0.71 


1.47 


2 


1.35787(2) 


0.364(1) 


0.376(1) 


3.116(8) 


0.89 


0.85 


1.37 


0.94 


3 


1.58951(2) 


0.377(1) 


0.386(1) 


3.212(8) 


0.57 


0.54 


1.07 


1.00 


/? = 0.20 





1.9141(1) 


0.0414(2) 


0.0884(2) 


1.458(9) 


0.55 


0.52 


0.47 


0.81 


1 


2.12263(9) 


0.0884(2) 


0.1106(3) 


2.14(1) 


0.41 


0.40 


1.38 


0.98 


2 


2.33075(8) 


0.1109(2) 


0.1178(3) 


2.40(1) 


0.43 


0.34 


1.02 


1.67 


3 


2.52196(8) 


0.1186(3) 


0.1150(2) 


2.48(1) 


0.56 


0.52 


1.36 


1.17 



TABLE I: Parameters of effective quantum mechanics that describes the data for the L = 2 lattice at 
/3 = 0.43. 

and Lt is sufficiently large. The corresponding effective parameters are tabulated in Tabs. |Ill |ni] 
andHYl The fits always give reasonable / DOF, which are shown in the last four columns, one 
for each observable. We note that as j3 becomes smaller, jif^ becomes larger while jii'^ — jj!f^ 
becomes smaller. This is the reason it becomes difficult to match the data to an effective quantum 
mechanics description at small (3 without going to very large Lt. Note also that the value of nf^ 
has approximately reached the thermodynamic limit dX 13 = 0.2 for L = 8. We plot our four 
observables near /ic"'' at L = 16 and /3 = 0.43 in Fig.[6]and at L = 6 and /3 = 0.2 in Fig.|7]along 
with the fits. 
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N 


L 




Po 


Pi 


Xo 


xVdof 





4 


1.85850(6) 


0.00100(3) 


0.02062(7) 


1.55(2) 


1.40 


1.50 


1.42 


0.85 





6 


1.85794(2) 


0.00002(1) 


0.00888(3) 


1.62(2) 


0.72 


1.05 


0.97 


0.98 





8 


1.85798(3) 


0.000001(1) 


0.00500(4) 


1.5(2) 


0.58 


0.56 


0.92 


0.39 





12 


1.85801(3) 


_ 


_ 


_ 


1.30 








1 


4 


1.90898(5) 


0.02018(8) 


0.03768(9) 


2.98(3) 


1.44 


1.61 


1.64 


1.41 


1 


6 


1.87771(2) 


0.00881(4) 


0.01731(6) 


3.04(6) 


1.07 


0.83 


0.56 


0.81 


1 


8 


1.86810(1) 


0.00497(3) 


0.00987(5) 


3.1(1) 


0.82 


0.77 


0.74 


1.32 


2 


4 


1.95988(4) 


0.0371(1) 


0.0515(1) 


4.05(3) 


1.51 


1.31 


1.62 


1.44 


2 


6 


1.89811(2) 


0.01732(6) 


0.02494(7) 


4.33(6) 


0.91 


0.82 


1.53 


0.48 


2 


8 


1.87866(1) 


0.00986(6) 


0.01446(4) 


4.4(2) 


0.61 


0.59 


0.63 


0.43 


3 


4 


2.01050(4) 


0.0512(1) 


0.0631(1) 


5.00(4) 


1.53 


1.46 


0.97 


1.40 


3 


6 


1.91911(1) 


0.02486(7) 


0.0323(1) 


5.65(7) 


1.59 


1.38 


0.92 


0.40 


3 


8 


1.88957(1) 


0.01451(5) 


0.01897(6) 


5.7(2) 


0.58 


0.55 


0.74 


0.42 



TABLE U: Effective Quantum Mechanics parameters near the various particle number transitions from N 
to AT + 1 at /3 = 0.20 and various values of the spatial size L. 
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N 


L 


P^ 


Po 


Pi 


Xo 


xVdof 





4 


0.45019(3) 


0.0733(6) 


0.1307(5) 


4.50(3) 


n 49 




Q9 


1 00 
1 .uu 





6 


0.35692(3) 


0.0330(2) 


0.0766(3) 


5.83(5) 








U.J J 


1 91 





8 


0.32270(3) 


0.0163(1) 


0.0489(2) 


6.3(1) 




1 01 


60 

U.UU 


49 





10 


0.30837(2) 


0.0081(2) 


0.0334(3) 


6.5(2) 


0.39 


0.43 


0.64 


0.99 





12 


0.30207(2) 


0.0043(3) 


0.0236(4) 


7.1(2) 


0.54 


0.52 


0.84 


0.84 





16 


0.29789(2) 


0.0013(4) 


0.0135(4) 


7.0(4) 


0.52 


0.52 


0.50 


0.48 





20 


0.29704(7) 


0.00039(5) 


0.0083(2) 


7.5(6) 


0.90 


0.80 


1.15 


0.36 


1 


4 


0.66047(2) 


0.1311(6) 


0.1670(6) 


6.07(3) 


1 04 


1 10 


66 

U.UU 


73 

U. / J 


1 


6 


0.48596(2) 


0.0763(3) 


0.1060(3) 


8.56(6) 




63 


1 93 


^7 

U.J / 


1 


8 


0.40856(2) 


0.0491(2) 


0.0728(2) 


10.2(1) 


0.38 


0.36 


1.13 


0.43 


1 


10 


0.36818(2) 


0.0336(2) 


0.0525(4) 


11.3(2) 


0.71 


0.84 


0.80 


1.55 


1 


12 


0.34516(1) 


0.0242(2) 


0.0387(2) 


11.9(2) 


0.59 


0.57 


0.90 


0.51 


1 


16 


0.32214(1) 


0.0136(2) 


0.0236(2) 


12.1(4) 


0.59 


0.60 


0.80 


0.90 


2 


4 


0.81586(2) 


0.1676(6) 


0.1944(6) 


7.10(3) 


"^4 




Q4 


1 00 
1 .uu 


2 


6 


0.58498(1) 


0.1057(3) 


0.1295(3) 


10.35(6) 


u.ou 


89 


78 
u. / o 


74 

U. / 


2 


8 


0.47730(1) 


0.0727(2) 


0.0922(2) 


12.9(1) 


0.65 


0.68 


0.65 


0.74 


2 


10 


0.41819(1) 


0.0522(3) 


0.0686(3) 


14.3(2) 


0.70 


0.66 


1.18 


0.76 


2 


12 


0.38263(1) 


0.0391(3) 


0.0525(3) 


16.1(3) 


0.78 


0.75 


1.06 


0.89 


2 


16 


0.34468(1) 


0.0233(2) 


0.0335(2) 


17.7(5) 


0.39 


0.37 


0.78 


0.75 


3 


4 


0.94318(2) 


0.1944(7) 


0.2150(5) 


7.94(3) 


50 


53 


99 


1 10 


3 


6 


0.66810(1) 


0.1293(3) 


0.1478(3) 


11.93(6) 


0.54 


0.87 


0.99 


0.48 


3 


8 


0.536227(9) 


0.0916(2) 


0.1079(3) 


15.1(1) 


0.47 


0.47 


1.24 


0.98 


3 


10 


0.46209(1) 


0.0679(3) 


0.0820(3) 


17.6(2) 


0.58 


0.57 


0.93 


0.59 


3 


12 


0.41633(1) 


0.0526(3) 


0.0637(4) 


19.2(3) 


0.46 


0.44 


0.80 


0.37 


3 


16 


0.365744(9) 


0.0328(3) 


0.0412(2) 


22.8(5) 


0.80 


0.76 


0.58 


1.08 



TABLE III: Effective Quantum Mechanics parameters near the various particle number transitions from N 
to A'^ + 1 at /3 = 0.43 and various values of the spatial size L. 
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N 


L 


/^^ 


Po 


Pi 


Xo 


xVdof 





4 


0.20651(2) 


0.1868(7) 


0.2173(7) 


7.01(3) 


fO 




n 79 


1 1 8 





6 


0.09166(2) 


0.1717(5) 


0.1847(5) 


13.39(9) 


xj.yjy 




1 91 







8 


0.049211(9) 


0.1691(5) 


0.1732(5) 


22.0(1) 


0.82 


0.79 


1.08 


0.80 





10 


0.030691(8) 


0.1677(7) 


0.1700(6) 


32.5(3) 


0.54 


0.51 


0.40 


0.64 





12 


0.02112(1) 


0.1677(4) 


0.1681(4) 


45.1(7) 


0.63 


0.60 


1.06 


1.10 





16 


0.011786(8) 


0.1670(4) 


0.1670(4) 


77.3(8) 


0.68 


0.63 


0.88 


0.74 


1 


4 


0.43892(2) 


0.2162(7) 


0.2440(7) 


7.76(3) 


u. o / 


n 8^ 


U. J>U 


n 68 


1 


6 


0.23271(2) 


0.1832(5) 


0.1997(5) 


14.00(9) 




n 4.0 


n 47 


1 1 9 


1 


8 


0.138584(9) 


0.1737(4) 


0.1801(5) 


22.3(1) 


0.72 


0.69 


1.23 


0.82 


1 


10 


0.090004(8) 


0.1699(6) 


0.1729(6) 


32.9(3) 


0.77 


0.74 


0.63 


1.5 


1 


12 


0.06268(1) 


0.1680(4) 


0.1704(5) 


45.6(7) 


1.50 


1.42 


0.52 


0.66 


1 


16 


0.035256(8) 


0.1672(4) 


0.1673(4) 


77.0(8) 


0.82 


0.78 


0.94 


0.70 


2 


4 


0.61012(2) 


0.2447(7) 


0.2638(7) 


8.40(3) 


u. ou 


n 8^ 


1 47 


1 0*^ 


2 


6 


0.34520(1) 


0.1986(9) 


0.2131(9) 


14.99(8) 


u. J>z. 




1 1 1 


84 


2 


8 


0.216371(9) 


0.1813(6) 


0.1891(6) 


23.1(1) 


0.79 


0.75 


0.79 


0.82 


2 


10 


0.145056(7) 


0.1729(6) 


0.1789(6) 


33.4(3) 


0.88 


0.83 


1.12 


0.94 


2 


12 


0.10266(1) 


0.1701(4) 


0.1734(4) 


44.8(7) 


0.72 


0.68 


0.49 


1.02 


2 


16 


0.058435(8) 


0.1680(4) 


0.1687(4) 


78.4(8) 


0.34 


0.32 


0.82 


0.85 


3 


4 


0.74886(1) 


0.2639(7) 


0.2824(7) 


8.98(3) 


n 79 


6Q 




46 


3 


6 


0.43967(1) 


0.2123(9) 


0.2280(9) 


15.74(8) 




n 40 


1 1 6 
1. ID 


n so 


3 


8 


0.284653(8) 


0.1900(6) 


0.1975(6) 


24.1(1) 


0.72 


0.69 


1.25 


0.79 


3 


10 


0.195544(7) 


0.1780(7) 


0.1838(6) 


34.8(3) 


0.63 


0.60 


1.09 


0.93 


3 


12 


0.14064(1) 


0.1729(4) 


0.1762(4) 


46.1(7) 


1.04 


0.98 


0.58 


1.50 


3 


16 


0.081185(8) 


0.1689(4) 


0.1698(4) 


78.3(8) 


0.88 


0.83 


0.66 


1.39 



TABLE rV: Effective Quantum Mechanics parameters near the various particle number transitions from N 
to + 1 at = 0.50 and various values of the spatial size L. 
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FIG. 3: The average particle number and the condensate susceptibility as a function of fi near the transition 
between = and = 1 at L = 2 and /? = 0.43 for different values of Lt. The solid lines show the fit 
of data to the effective quantum mechanics description discussed in the text. 
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FIG. 4: The average particle number as a function of for at /? = 0.43 (left) /3 = 0.50 (right) for different 

values of L. When /3 = 0.43 the data shown is for = 100 at L = 2 8, Lt = 200 at L = 12 and 

Lt = 300 at L = 16. When /3 = 0.50 the data shown is for = 100 at L = 2 6, Lt = 160 at L = 8, 

Lt = 200 at L = 12 and Lt = 300 at L = 16. 
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FIG. 5: The average particle number as a function of ^ near the transitions between and + 1 for 
iV = 1, 2, 3 at L = 8 and /3 = 0.43 (top), (3 = 0.5 (center) and (3 = 0.2 (bottom).The solid hnes are fits to 
the data. 
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FIG. 6: The four observables as a function of /x near the transition between = and = 1 at L = 16 
and /3 = 0.43 at two different values of Lt. The solid lines show the fit of data to the effective quantum 
mechanics description discussed in the text. 
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FIG. 7: The four observables as a function of /i near the transition between = and A'^ = 1 at L = 6 and 
/3 = 0.2 for Lt = 540 and 1080. The solid lines show the fit of data to the effective quantum mechanics 
description discussed in the text. 
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VI. THERMODYNAMIC LIMIT 



Using the results of the previous section it is tempting to extrapolate to the thermodynamic 
limit. However, in order to accomplish this task it is important to know how the effective quantum 
mechanical parameters depend on L. This dependence is non-universal in general and close to 
a critical point will depend on the nature of the phase transition. Assuming the phase transition 
is second order, close to the critical chemical potential where the density can be made arbitrarily 
small, we expect universal features to emerge. For example, when the particles have a purely 
repulsive interaction, the ground state energy of particles is always less than the corresponding 



energy of + 1 particles 13 ID . Based on the results of the previous section this scenario seems to 
be valid in the current model. Indeed the particle number always increases by one as we increase 
ji at every fixed value of L. The superfluid density ps also behaves like p. Thus, we conclude that 
at /i = pi^ in the thermodynamic limit, there is a second order transition to a superfluid phase. 
Based on this, below we discuss the extrapolations to the thermodynamic limit. 

First we consider /3 = 0.43 where the low energy physics contains massive bosons with repul- 
sive interactions. Then, the quantity pf"* is simply the mass of the particle M{L) at a finite L. This 
mass can be obtained by fitting the the temporal two-point correlation function 

G{t) = J2 (e*^^^e-*^^^y (22) 

computed at /i = 0, to the form G{t) ~ exp(— M(L)t) for values of t ^ Lt/2. In the definition of 
G{t), y± and x± represent lattice sites at temporal slices and t respectively. We have computed 
M{L) using this method and indeed we find excellent agreement with pl^^ at all the three values 
of /3. This means the true mass of the particle must be 

M = lim p^°\ (23) 

We can reverse this argument and obtain /ic°^ in the thermodynamic limit by simply measuring the 
mass of the particle at p = 0. Of course this result is not general and is valid only in the present 
study where there is clear evidence that the particles re pel each other. In order to extract M in the 



massive phase (/3 < (3c) we can use Liischer's formula [|32|l extended to two spatial dimensions. 



/i° ^ M + MiQ-^^. (24) 

At /3 = 0.43 we find that fits weU to this form and gives M = 0.29680(4), Mi = 0.59(2) and 
m = 0.393(4) with a small x^/DOF. The data and the fit are shown in the left plot of Fig. [8l 
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FIG. 8: The finite size scaling of tlie N-particIe energy levels as a function of the spatial lattice size L at 
/? = 0.43 (left, massive phase) and j3 = 0.50 (right, superfluid phase). 

The spatial size dependence of the energy of N particles in three spatial dimensions has been 
calculated using models of quantum mechanics JSJ, Recently, this dependence was also 
computed using effective field theory [Ssj]. In the special case of two particles the problem was 
also solved in a general massive quantum field theory in three spatial dimensions [|36|] . All these 



studies indicate that the ground state energy of N particles satisfies the relation e'^'^ — e'^'' oc 
A^(A^ + 1)/L^- For > 1, remember that /ic^^ is the difference in the ground state energies 
of + 1 particles and N particles. Extending the known results to two spatial dimensions and 
making the assumption that the particle density in the thermodynamic limit close to the critical 
point is of the form p ~ c(/i — /ic°^) at leading order in the superfluid phase, we expect 

(iV+1) 



cL2 



(25) 



for sufficiently large L and A^. Figure [8] shows that our data is described reasonably well by this 
equation. In the left plot of Fig. [8] we show the values of fi[^\L) obtained from table Hill The 
solid lines show the dependence of fii^\L) on L as described by eq. (l25l) with c ~ 0.18, 0.16 and 
0.15 at = 1, 2 and 3 respectively. Clearly, for large values of L the solid lines pass through the 
data. The value of c changes slightly since A^ is small. Unfortunately a fit of our data to eq. (1251) 
yields a large x^/DOF. We believe this is due to the fact that our data has very small errors and 
hence is sensitive to higher order corrections which we do not know analytically at the moment in 
two spatial dimensions. 
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When /3 = 0.5 we are in the superfluid phase and the U{1) particle number symmetry is 
spontaneously broken. One then expects the low energy spectrum at finite volumes to be governed 
by 0(2) chiral perturbation theory. Based on this we again expect jj,i^\L) to be described by 
eq. (|25l) but with = 0. While our data is again consistent with these expectations (see right 
plot of Fig. [8]), without keeping higher order 1/L corrections the fits again give a large x^/ DOF. 
The solid lines in Fig. [8] describe eq. ^ with c = 0.33, 0.22, 0.2, 0.2 for = 0, 1, 2, 3. 



VII. PHASE DIAGRAM 

The phase diagram of the 0(2) non-linear sigma model is an interesting research topic in itself. 
While the complete phase diagram requires more work, our results above allow us to compute the 
location of the transition line between the normal phase and the superfluid phase. In particular 
the value of /ic°^ as a function of /3 determines this line. Based on the evidence at /3 = 0.43 and 
0.20 we predict that = M for all values of /3 < /3c- The coordinates of the transition line 
are tabulated in Tab. |V]and shown on the phase diagram in Fig. [9l We expect this transition to 
be second order in the mean field universality class with logarithmic corrections except at yU = 
where it is governed by the 3(i XY universality class. Thus, when /3 < f3c and (/3c — l3)/l3c <^ 1 
we must have /i° oc [(/3c - P)/Pc]'' where u ^ 0.671 [37]. 
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0.43 


0.40 


0.35 


0.30 


0.25 


0.2 


0.15 


0.10 


0.07 


0.04 




0.29678(3) 


0.530(4) 


0.859(2) 


1.172(7) 


1.505(9) 


1.85801(3) 


2.267(3) 


2.783(3) 


3.210(3) 


3.829(4) 



TABLE V: The values of n^c^ obtained by assuming that it is equal to the L ^ oo limit of M{L) at /i = 
as discussed in the text. 



In principle there could be other interesting phases at larger values of ji which we cannot rule 
out based on the current work. Since we have seen the particles have a repulsive interaction, an 
interesting possibility is the existence of a solid phase or a super-solid phase 113 811 . However, there 



4011 and we do not know if these rule out such 



are stringent constraints for super-solids to arise [I39I .I' 
a phase in the current model. In any case if a transition to a solid phase exists, it will most likely be 
first order similar to the solid-liquid phase transitions in materials and will occur at densities where 
the lattice structure may become important. These transitions can also be studied efficiently with 
the worm algorithm. We postpone these studies for the future, but have speculated the possibility 
of a solid phase in Fig. |9l 
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p 

FIG. 9: The phase diagram in the /3 vs. fi plane. The circles show the value of fi^^ as a function of /? given 
in table |Vl The solid line that connects these points forms the phase boundary between the normal phase 
and the superfluid phase. This transition is second order. Since the particles repel each other we speculate 
that at higher densities a first order transition (dashed Une) may separate the superfluid phase from a solid 
or a super-solid phase. 
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Appendix A: The Worm Algorithm 



The worm algorithm for the partition function described by eq. |4]can easily be constructed 
using ideas from II101I141 I20I1. Here we outline the essential steps of the update for completeness. 
Each worm update is as follows 

1. We pick a random point x on the lattice. We will also call this site xgrst- We set a counter 
c = 0. 

2. We pick at random one of 2d neighbors x + ce, a = ±1, ±2, ±d of the site x. 

3. Let k be the current on the bond connecting x and x+a. If a is positive then with probability 

we change A; to A; + 1 and move to the neighboring site x + a. If a is negative then with 
probability 

/fc-i(^)e-^'^".« 

we change kto k — 1 and move to the neighboring site x + a. Otherwise we stay at site x. 

4. We set c = c + 1. If a; = Xfirst we stop and complete one worm update. Otherwise we go to 
step 2 and repeat the above steps. 

It turns out that x is given by the average of c after many worm updates. The other observables 
are measured on each world-line configuration and averaged over the ensemble generated by the 
worm algorithm. 



Appendix B: Tests of the Algorithm 

We have verified our algorithm and code by both solving the model exactly on a 2 x 2 lattice as 
well as comparing with the available results in the literature for /i = in three dimensions. In this 
section we describe some of these tests. First, we compare the results of the various observables 
computed using the directed path algorithm with the exact results on a 2 x 2 lattice. The comparison 
is shown in Table IVTl Since space and time are symmetric we expect = k at /i = 0. Our results 
reflect this fact. 
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/3 


X 




P 




K 




Ps 




fi = 0.0 


0.1 


1.2207 


1.2206(1) 





2(2) X 10-^ 


0.01005 


0.01001(4) 


0.01005 


0.01008(4) 


0.2 


1.4831 


1.4830(2) 





-5(4) X 10-5 


0.04064 


0.04063(8) 


0.04064 


0.04076(8) 


0.5 


2.3838 


2.3835(4) 





1(9) X IQ-^ 


0.2526 


0.2527(2) 


0.2526 


0.2526(2) 


1.0 


3.2730 


3.2733(4) 





2(2) X 10-^ 


0.7796 


0.7796(4) 


0.7796 


0.7794(4) 


5.0 


3.8728 


3.8733(3) 





1(1) X 10-^ 


4.809 


4.806(4) 


4.809 


4.813(4) 


fi = 0.5 


0.1 


1.2362 


1.2361(1) 


0.00590 


0.00592(2) 


0.01559 


0.01560(5) 


0.01007 


0.01013(4) 


0.2 


1.5188 


1.5188(2) 


0.02374 


0.02378(5) 


0.0640 


0.0640(1) 


0.04100 


0.04105(7) 


0.5 


2.4692 


2.4693(4) 


0.1429 


0.1430(1) 


0.4341 


0.4344(3) 


0.2581 


0.2581(2) 


1.0 


3.3245 


3.3247(4) 


0.4190 


0.4193(2) 


1.6553 


1.6564(10) 


0.7867 


0.7861(4) 


II = 1.0 


0.1 


1.2861 


1.2861(2) 


0.01809 


0.01812(4) 


0.03851 


0.03857(8) 


0.01018 


0.01024(4) 


0.2 


1.6316 


1.6317(2) 


0.07160 


0.07162(8) 


0.1637 


0.1638(2) 


0.04241 


0.04242(8) 


0.5 


2.6899 


2.6900(4) 


0.3865 


0.3868(2) 


1.2530 


1.2539(8) 


0.2747 


0.2748(2) 


1.0 


3.4317 


3.4317(4) 


1.0053 


1.0055(3) 


5.487 


5.488(3) 


0.8032 


0.8033(4) 



TABLE VI: Checks of the observables with exact solution on 2 x 2 lattices. 



Extending the code from two dimensions to three dimensions is trivial and the chance for 
mistakes is rather small. However, we have tested the code at least at yu = using the results 
from previous work. Here we compare results for x obtained using the worm algorithm with that 
obtained using the microcanonical improved Metropolis (MM) algorithm and the available results 
in the literature using the cluster method [28]. The comparison is shown in Table lVIII The reason 
for us to choose (3 = 0.45421 is because this is known to be the critical value of the coupling 
where the theory undergoes a phase transition from a normal phase to a superfluid phase. At the 
critical coupling we expect x L"'^'^. A fit of our data to this form yields the value of 7/z/ = 1.99 
as expected from n2m . and is shown in figure [TOl 
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Susceptibility per unit site using various metliods 
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FIG. 10: Plot of X vs L at /3 = (3^. 





L 


Worm 


MM 


Cluster 


0.45421 


4 


19.17(3) 


19.24(13) 


19.15(5) 


0.45421 


8 


77.8(2) 


76.9(5) 


77.9(3) 


0.45421 


16 


310(1) 


313(2) 


313(2) 


0.45421 


32 


1221(18) 


1228(7) 


1226(13) 


0.45421 


48 


2713(67) 


2750(27) 


2719(68) 


0.01 


8 


1.0304(2) 


1.03(46) 




0.1 


8 


1.3976(9) 


1.40(6) 




1.0 


8 


387.2(3) 


387.12(3) 





TABLE VII: Comparison of the condensate susceptibility x with results from the worm algorithm, the 
Metropolis+Microcanonical(MM) update and Wolff Cluster update on lattices at different values of L 
and p. 
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Appendix C: Raw Data for comparison at /x 7^ 



Here we give some raw data obtained using our algorithm at non-zero chemical potential for 
comparison with other methods like the complex Langevin method which are being developed to 
solve the sign problem present in the conventional formulation. These values can also serve as a 
check for future work. 



P 


L 


Li 


/' 


P 


/,■ 


Ps 


\ 


0.5 


12 


144 


0.035 


0.0625(1) 


0.932(2) 


0.1685(7) 


3960(11) 


0.43 


4 


100 


0.446 


0.0248(4) 


2.48(4) 


0.0972(8) 


239(2) 


0.43 


16 


960 


0.297 


0.00120(3) 


1.15(3) 


0.0050(1) 


3422(31) 


0.30 


8 


8 


1.20 


0.0450(3) 


1.39(1) 


0.0099(1) 


34.8(2) 


0.30 


8 


64 


1.20 


0.0246(2) 


3.17(3) 


0.0164(1) 


250(2) 


0.20 


6 


180 


1.855 


0.0104(2) 


1.94(4) 


0.00326(6) 


188(2) 



TABLE VIII: Raw data obtained using the worm algorithm at random values of the parameters 



Lt 




P 


K 


Ps 


X 


320 


0.404 


0.01850(8) 


7.76(8) 


0.0535(4) 


1509(13) 


320 


0.406 


0.02038(9) 


9.57(9) 


0.0562(4) 


1661(12) 


320 


0.408 


0.0227(1) 


11.8(1) 


0.0597(4) 


1752(15) 


320 


0.41 


0.0252(1) 


14.2(1) 


0.0646(4) 


1766(13) 


320 


0.412 


0.02745(8) 


16.35(8) 


0.0668(5) 


1658(12) 


320 


0.414 


0.02890(6) 


17.75(6) 


0.0692(4) 


1503(11) 



TABLE IX: Raw data at /3 = 0.43, L = 8 near the transition from N=l to N=2. 
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